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TRAN SVERSE VIBRATIONS OF HOLLOW THIN-WALLED 
CYLINDRICAL BEAMS! 


By BERNARD BuUDIANSKY and Epwin T. Krouszpwsxki 


SUMMARY 


The variational principle, differential equations, and bound- 
ary conditions considered appropriate to the analysis of trans- 
verse vibrations of hollow thin-walled cylindrical beams are 
shown. General solutions for the modes and frequencies of 
cantilever and free-free cylindrical beams of arbitrary cross sec- 
tion but of uniform thickness are given. The combined influence 
of the secondary effects of transverse shear deformation, shear lag, 
and longitudinal inertia is shown in the form of curves for 
cylinders of rectangular cross section and uniform thickness. 
The contribution of each of the secondary effects to the total 
reduction in the actual frequency is also indicated. 


INTRODUCTION 


The elementary theory of bending vibration is often in- 
adequate for the accurate calculation of natural modes and 
frequencies of hollow, thin-walled cylindrical beams. Such 
secondary effects as transverse shear deformation, shear lag, 
and longitudinal inertia, which are not considered in the 
elementary theory of lateral oscillations, can have appre- 
ciable influence, particularly on the higher modes and 
frequencies of vibration. The effects of transverse shear 
deformation and of rotary (rather than longitudinal) inertia 
have been studied by many on the basis of the original inves- 
tigations of Rayleigh (ref. 1) and Timoshenko (ref. 2). 
Anderson and Houbolt (ref. 3) have presented a procedure 
for including the effects of shear lag m the numerical calcu- 
lation of modes and frequencies of box beams of rectangular 
cross section. However, there does not appear to exist a 
general solution for the vibration of hollow beams that in- 
corporates the influence of all the secondary effects men- 
tioned. 

The purpose of the present report is threefold: First, to 
exhibit the variational principle, differential equations, and 
boundary conditions appropriate for the analysis of the 


uncoupled bending vibration of hollow thin-walled cylindrical | 


beams; second, to give general solutions for cantilever and 
free-free cylinders of arbitrary cross section but of uniform 
thickness; and finally, to show quantitatively the influence 


of the secondary effects by means of numerical results for 
hollow beams of rectangular cross section of various lengths, 
widths, and depths. 


SYMBOLS 
A cross-sectional area 
Ax Fourier coefficient 
Ag effective shear-carrying area - 
B; parameter defined in equation (30) 
C constant 
B- modulus of elasticity 
G shear modulus of elasticity 
I moment of inertia . 
K geometrical parameter defined in equation. (29) 
L length of cantilever beam, half-length of free- 
free beam 
Ni parameter defined in equation (38) 
T maximum kinetic energy 
U maximum strain energy 
a half-depth of rectangular beam 
b half-width of rectangular beam 
Amn, On Fourier series coeficients 
i,j, m,n integers 
; [uD 
kp frequency coefficient, w mt 
A EI 
8 coefficient of shear rigidity, 7 A.G 
ny coefficient of rotary inertia, Liz 
p. perimeter of cross section 
8 distance along Penpnery of cross section (see 
fig. 1) 
t wall thickness 
u(z,8) longitudinal displacement in z-direction 
w(x) vertical displacement in y-direction 
2 longitudinal coordinate 
y ` - vertical coordinate ° i 
y y-coordinate of center of gravity of cross section 
Ves shear strain 


1 Supersedes NAOA TN 2682, “Transverse Vibrations of Hollow Thin-Walled Cylindrical Beams” by Bernard Budiansky and Edwin T. Kruszewski, 1952, 
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€z longitudinal strain 

0 inclination of normal with vertical (see fig. 1) 

A Lagrangian multiplier 

m mass of beam per unit length 

p mass density of beam 

o longitudinal direct stress 

T shear stress 

w natural frequency of beam 

Wo natural frequency of beam calculated from ele- 
mentary beam theory 

Say Kronecker delta (1 if <= 9; 0 if 757) 

p constraining relationship 


BASIC EQUATIONS 


Assumptions.—The problem to be considered is that of 
the natural bending vibration of a thin-walled hollow 
cylindrical beam whose cross section is symmetrical about at 
least one axis (see fig. 1). The transverse vibration is sup- 
posed to take place in the direction of this axis of symmetry 
of the cross section so that no torsional oscillations are 
induced. 

In the present analysis, the following simplifications are 
introduced: 

(a) Changes in the size and shape of the cross section are 
neglected. yo 

(b) Stress and strain are assumed to be uniform across the 
wall thickness. ` 

(c) The small effect of circumferential stress upon longi- 
tudinal strain is neglected. ` 

In accordance with statements (a) and (b), the distortions 
of the vibrating beam are completely described by the 
vertical displacement w(x) of a cross section and the longi- 


y 
| V6.9. 
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a) Typical cross section. Sign conventions. 
(3) Cantilever beam. ts Free-free beam. 


Figure 1.—Coordinate systems and sign conventions. 
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tudinal displacement w(z,s) of each point of the median line 
of the beam wall. ` 


The longitudinal-and shear strains are given in terms of 
u(z,s) and w(z) as 


Ow 
rs 1) 
and 
ou, dw . 
Ya=5, tq, #0 ð (2) 


g=E o l (3) 
and 
òu , dw . 
Tu= G ae Er sın o) (4) 


where 0 is the inclination of the normal with the vertical (see 
fg. 1). i n 
In elementary beam theory, where the effects of all shear 


- distortion are neglected, the longitudinal distortion u(z,8) is 


related to the vertical displacement w(x) by 
d 
u (£,8)=Y—Y) o 


where ¥ is the y-coordinate of the center of gravity of the 
cross section. In the present report, however, w(z,s) is 
allowed to be perfectly general, so that shear distortions (and 
consequently the so-called shear-lag and transverse-shear- 
deformation effects) are fully taken into account. Further- 
more, because cross sections are not constrained to remain 
plane, the inertia effect associated with motion in the 
longitudinal direction is more properly designated as the 
effect of longitudinal inertia than the effect of rotary inertia. 

' Variational principle and geometrical boundary condi- 
tions.—The variational equation to be written is appropriate 
to beams whose ends are either fixed, simply supported, or 
free. For some such beam vibrating in a natural mode, the 
maximum strain energy is, 


1 (F? our? 1 (2 ou, dw. 2 
U=} Í, PE( Se) tdede+5 Í ha t Fein 0) dedz 
(5) 


where w(x,8) and w(x) are the amplitudes of displacement for 
the particular mode considered. The maximum kinetic 
energy is ' 


1 fz a 1 (2 
=i f O ptutwtds dett f pturutds de (0) 


where w is the natural frequency of the mode under considera- 
tion and p is the mass density of the beam. The second 
term in equation (6) constitutes the contribution of longi- 
tudinal inertia to the kinetic energy. 

A natural mode of vibration must satisfy the variational 
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equation 


ô(U —L)=0 (7) 


where the variation is taken independently with respect to 
u(z,s) and w(x) and with the provision that both u(z,s) and 
w(x) must satisfy the geometrical boundary conditions of 
the problem; furthermore, w(z,s) must be periodic in the 
coordinate 3 with a period equal to the perimeter p. The 
geometrical boundary conditions are w=0 and u=0 at a 
fixed end and only w=0 at a simply supported end. Ata 
free end no geometrical boundary conditions are imposed. 
Differential equations and natural boundary conditions.— 
Equations (5), (6), and (7) in conjunction with the usual 
procedure of the calculus of variations yield the following 
simultaneous integrodifferential equations for u and w: 


Ti u, Ò Ow, dw. 
it Se Ds al? Os Je de dz" p) |+ptat “aa (8) 
p Gt Kaa = sin + Te gin? 8 0) ds+ po*w=0 (9) 


where 


(10) 


w= pt ds 


and the boundary equations at each end of the beam are 


GrH) mde=o aw 


oe U1? sin 0) sin 0 ds ow—0 


fat Gt 

At a fixed end, both boundary equations (11) and (12) are 

satisfied by virtue of the fact that the geometrical boundary 

conditions require that both éu and ôw be zero. At a simply 

supported end éw=0, but, since du(z,s) is perfectly arbitrary, 
the variational process forces the equality 


(12) 


OU 
Et F. 


(13) 
Finally, at a free end, since there are no geometrical con- 
straints, both ôu and ôw are arbitrary and hence the varia- 
tional process forces, in addition to equation (13), the 


equality 
dw . ; 
han os ae sin ) i sin 0 ds=0 


Equations (13) and (14) constitute so-called “natural 
boundary conditions” because they are automatically satis- 
fied as the result of a variational process. Equation (13) is 
recognized as the condition of zero longitudinal direct stress 
while equation (14) simply stipulates that the total vertical 
shear force vanish. 


Thus to summarize, the appropriate boundary conditions 
required for the solution of equations (8) and (9) are 


(14) 
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Fixed end: 
" w=0 
u=0Q0 
Simply supported end: 
aw=0 
qn, OU 7 
Et = 0 
Free end: 
ou , dw . 


pals 


TaT Ey sin a) sin 0 ds=0 


OU 
Et pO 


The integrodifferential equations (8) and (9), which specify 
equilibrium in the longitudinal and transverse directions 
respectively, can, of course, be written directly without 
recourse to the variational principle. 


GENERAL SOLUTIONS FOR CYLINDERS OF UNIFORM WALL 
THICKNESS 


The following exact solutions for cylinders of uniform wall 
thickness are carried out by means of Fourier series in con- 
junction with the application of the variational condition 
(eq. (7)). This procedure, which does not require explicit 
consideration of the natural boundary conditions, was 
believed to be more expedient than a direct attack upon the 
simultaneous integrodifferential equations (8) and (9) and 
all their associated boundary conditions. 

Cantilever beam.—The geometrical boundary conditions, 
for a cantilever beam, as previously shown, are 


w(0)=u(0,8)=0 

(see fig. 1). Appropriate assumptions for the pean l 
w(x) and u(z,8) are 

w(x)=C+ 2. ba COS a (15) 


and 
2NTE 


ulx,8)= >> a Amn sin — 


1 
m=1,3,5 n=0,1,2 2L © (16) 


The condition u(0,s)=0 is satisfied by each term of equa- 
tion (16); the condition 


wO)=C+ > b,=0 (17) 


is introduced into the variational procedure by means of the 
Lagrangian multiplier method. The choice of the particular 
trigonometric functions used in the Fourier series (15) and 
(16) was guided by consideration of the orthogonality 
required for the simplification of expressions in the strain 
energy. The constant C is needed in the expression for 
w(x) in order that w(Z) be unrestricted. 


ee ee ee ee a e -e ee va m - e a 
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Using equations (15) and (16) in equations (5) and 8) 
yields 
U—T= 


= mar Marz Inrs? 


1 L œ 
z, prl, 2 sad aia“ 2L OL ae P neat 
| fi = = 2ne . Max . NTE 
a fa (È. nia P iia 2L i P T 
i æ NT . NTEN? 
Sn Ad BLS OE ds dx— 

L rad 
al rotor ( 2s ba cos rn ds dx— 

: Maz a san) 

=i, Ge (33, n=0,1,2 mE gT ao Sa a 


To make equation (18) stationary and at the same time 
satisfy the constraining relationship 


, p=0+ 3 b,=0 (19) 
nal,3,5 


it is sufficient to set 
SU —T— ro) =0 (20) 


where the variation is with respect to the a’s, 6’s, and C 
considered as independent variables; here ` is a Lagrangian 
multiplier. This variational process results in the following 
equations: 


ann’ A, = 


OUT) _, Oe ina? “E = 
Ob; Ob, P D mta i a 
po” Ei Te 22 pat 
=0 (i=1, 3,5, ...) (21) 
DU — br NL 27% \? i 
er sr) Ft doa t At (==) M aut 
qitim Åjbi a 2a +o = 0 
(@=1, 3,5, .. .) 
=0,1,2,...) (22) 
es D oo 
oC Page -ut Sy 22 (ay 2 D n— po LO—d 
s i (23) 
where 
2 : . 2NTSE 
A=2h sip. 0 sin. ds (24) 
P P - 


As=h t sin? 6 ds (25) 
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With the use of the nondimensional parameters 


pT 
k= ir o (26) 
EI 
k= 10E (27) 
I I 
; ke AD (28) 
167 
K=7 (29) 
and . 
; 2\3 
Bea? —keptheg? (2) (30) 
equations (21), (22), and (23) may be reduced to 
2 inr? L l /ir\? I 1 
Sela de dg tata (3) pa beg bake 
yF 2 TH En oD =0 (i=1,3,5,...) (81) 
(k BEHEA) (1 +803) diy +E fr Ajijb;=0 
(1==1,3,5,...) 
(7=0,1,2,...) (32) 
t= =] 2 i 
ket 21, em | Si Tb +hF0+2=0 (33) 
For 7=0, equation (82) becomes 
kg? | hatha? 5] Aip==0 Q=1, 3, 5, ea ) (34) 


Equation (34) is not coupled to any of equations (31) to 
(33). A given value of a corresponds to the amplitude of 
the ith mode of longitudinal oscillation, and if this value of 
Go is not equal to 0, then equation (34) simply gives the 
frequency of this longitudinal mode. Consequently those 
equations in equation (32) for values of 7=0 are not asso- 
ciated with transverse bending and so are ignored henceforth. 
For the remaining values of 7 (that is, 70) equation (32) 


- yields 
n E ore 
am Rae Ani ,,  @=1,3,5,...) 
a ks B+E? i (j=1,2,3,...) (35) 


Substituting the expression for a in equation (35) into 
equation (31) and solving for b; gives 


CDe ra 


b= N, 


' @=1,3,5,...) (86) 
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where 
ate age: Sar 
nate N16 bs? As * pci 143 (37) ` 
© 8 kè ke Be+Kn? 2P 
n=l,2,3 
In the appendix: this expression for N; is shown to be ‘ee to 
Ur p, Ur A 2 n 2 
Ni=3q Pi — Te As ks? TEA BEFEA -zts | (38) 


Since the series in equation (38) is considerably more quickly convergent than that in equation (37), equation (38) should be 


used in actual numerical calculations of N3. 


Substitution of equation (36) into equation (33) and the constraining-relationship equation (19) gives the following two 


homogeneous equations in C and 2: 


be [he A (2) we] iA 


E that 3 ( yr Zz na SAG HORN, Na 


z1 3 172° 


Finally the condition for a nontrivial solution for C and à gives the frequency equation 


bot [at 32 Ge) N 


n-i 2 1 
Lhe D, (—1) $ A 
which the frequency parameter kg must satisfy. Since the 
terms of the infinite series which appear in the frequency 
equation contain / itself, the roots of equation (40) are most 
conveniently found by trial. Fortunately the infinite series 
in equation (40) as well as the series in the definition of N; 
converge rapidly so that only a few terms are needed to 
evaluate them with sufficient accuracy. 


Once kg has been determined for a particular mode, the 
corresponding mode shape can be found by letting C=1 and 
solving cither of equations (39) for ` and then finally evaluat- 
ing b; and a successively from equations (36) and (35). 


Free-free beam—symmetrical modes.—lIf the origin of a 
free-free beam of length 2Z is taken at the midspan (see fig. 
1), the form of the Fourier series assumed for w(z) and u(z,8) 
when the beam is undergoing a symmetrical mode of vibra- 
tion may be exactly the same as that assumed for the 
cantilever beam of length Z (see eqs. (15) and (16)). The 
only difference in the ensuing calculations is that the con- 
straining condition (19) is not introduced. Consequently, it 


can be readily seen that the frequency equation for the 


39 
24 (— 1) * nt Na EI ae ( 2) 
(39b) 
tet B ODP 2 ot 
: nr Na 
= 0 (40) 
(e) iL. . 
na13,5 Na 


symmetrically -vibrating free-free beam is obtained from 
equation (39a) by setting `A=0 and is 


2 2~ 2\F 1 7_ 
Š ks Ea Pa (= y|=° 
After a particular root kg is found from equation (41), the 
shape of the corresponding symmetrical free-free mode may 
be obtained from equations (36) (with A=0) and equations 
(35). 

Free-free beam—antisymmetrical modes.—Consider a 
free-free beam of length 2C undergoing antisymmetrical 
vibrations. Explicit consideration need be given only to the 
right half of the beam (see fig. 1), and for this half-beam the 
only geometrical boundary condition that must be imposed 
is that w(0)=0. The spanwise displacement u(0,s) is unre- 
strained by virtue of antisymmetry. 


Appropriate assumptions for the displacements w(e) and 
uw(2,8) are then 


(41) 


w(x)= >) b, a 


(42) 
nm, 4,6 
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u(z,3)= >) = y ann cos oT cos = (4.3) 


m=(Q,2,4 n=l, 


The linear portion Cx of the expression for w(x) is needed in order to give the beam sufficient freedom at the tip (=L). 
The choice of the particular trigonometric function in the series expansion for u(z,3) was, as in the case of the cantilever 
beam, guided by consideration of the orthogonality required for the simplification of the expressions in the strain energy. 
The zeroeth term in the series for u(z,8) in the s-direction was omitted because it only leads to the frequency eo for 


longitudinal oscillations. 


Using equations (42) and (43) in equations (5) and (6) yields 


ST 


w Mr i 
a at, se To is CeeS OL. 


sin ô 
nw2,4,6 


| t +( > 5a ma COS LEF cos 
= a ee 
2Jo j mo024n=123 7 2L 


The variation of equation (44) with respect to the a’s, 
b’s, and C gives, after suitable simplification, 


i 


Aji7b,s=0 Q= 2, 4, 6, -) 
(= 1,2,3, se) (45) 


(Biks?+-K27%) a4j—Fe& aL 


| BrP kstkatks (2) |u £ Ayjo=0 


(j=1,2,3,...) 46) 


ks? D = (— 1254+ = ke CL= 0 (48) 


nm2,4,6 NE 


From equation (45) 


il al ) 49) 
Qisj=-De2) P22 Vt = 2, 4, 6, . 49 
l Bèk’ +E F . (j=1,2,3,.--) 


which, except for sign, is the same expression as that ob- 
tained for the cantilever and symmetrically vibrating free- 
free beams (eq. (35)). From equation (46) 


BE Ag 


t= gah FEP C (j=1,2,3,.. 


.) (50) 
Substitution of equation (49) into equation (47) gives 


~~ pn 2 OL GH2,4.6 
t ( ) ix N: (2 3 EI Yj» e 


Inre? Jedi 
Pp 


) 61) 


_ 2nes Nnr mre. Br 
isda an fal $ ess 2nrs 
ze) += Pay pa —d rs cos —— OT, sin 5 -} 


5a bo aT cos weto) | naif $ pta? T> 24, b,, Sin m2 +Ox) ds dz— 


(44) 


where N; is defined in equation (37). 

Substitution of equations (50) and (51) into equation (48) 
and simplification gives as the frequency equation for the 
antisymmetrically TOES HRE free-free beam 


2 kg i co A 2 1 
2 B ar? a Ti 
j 2s on) 8 Netw ae! a aha Boke Ra 3] : 


(52) 


After a particular value of kg is found from equation (52), 
the shape of the corresponding antisymmetrical free-free 
mode may be obtained by giving C the arbitrary value of 
unity and calculating the b's and a’s successively from 
equations (51), (50), and (49). 

Discussion of parameters.—The parameters entering in 
the frequency equations merit discussion. The unknown 
natural frequency is contained only in the A co- 


efficient ks, which is defined by the formula w= kp a 


and is in common use in beam-vibration analysis. The 
parameters ks and kp; are identical with the shear and 
inertia parameters defined in reference 4, which considers 
the effect of only transverse shear and rotary inertia on 
beam vibrations. The quantity As which appears in the 
present definition of kg is actually the effective shear-carrying 
area, when plane sections are constrained to remain plane; 
that is, when shear lag is neglected. The remaining param- 
eters appearing in the present derivation, namely, A/Ag, 
K, and A;, Ag, . . . are essentially shape parameters which 
actually depend only on the contour of the cross section; 
as shown in the appendix, 


As_1 < 2 
A 3,253" 

and 
wee SAS 


T Ag n5i23 n’ 


and the A,’s are simply the Fourier coefficients of the 
function sin 0, which is dependent only on the shape of the 
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cross section. These shape parameters are related to shear- 
lag effects and their interaction with transverse shear 
and longitudinal inertia. 

The effect of longitudinal inertia is associated with the 
parameter kpr. If the effect of longitudinal inertia is to be 
neglected, it is sufficient to set ‘rr equal to zero in the final 
frequency equation. If kp; is equal to zero, B; becomes in- 
dependent of ks. Appreciable simplification in a trial-and- 
error solution for the natural frequency then results since, 


with B; independent of kg, the infinite summation contained - 


in N; is also independent of ks and need be calculated only 
once for any particular beam. As is shown in the following 
section, the effect of disregarding the influence of longitudinal 
inertia may often be negligible. 

Without presentation of details, it may be mentioned that 
for the case of a circular cylinder, which has no shear lag, 
all the A,’s except A; vanish and the frequency equations 
(40), (41), and (52) may be put into closed forms identical 
to those given in reference 4. Again, if in the general fre- 
quency equations ks is set equal to zero, the equations may 
be put into closed forms equivalent to those of reference 4 
where only rotary inertia is considered. 


‘ 
t 


RESULTS FOR CYLINDRICAL BEAMS OF RECTANGULAR 
CROSS SECTION 


In order to show quantitatively the effects of shear lag, 
transverse shear deformation, and longitudinal inertia on 
the natural frequencies of hollow thin-walled cylindrical 
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beams, numerical calculations have been performed for 
cylinders of rectangular cross section oscillating as free-free 
beams. The calculations have been limited to symmetrical 
modes of vibration, and consequently the frequency equation 
(41) is applicable. For rectangular cross sections the quan- 
tity N: may be put into closed form as shown in the appendix, 
and this closed-form version of N; was used in the calcula- 
tions. A value of E/G equal to 2.65 appropriate for 
aluminum alloys) was assumed. 

The results of these calculations are shown in izis 2-3, 
and 4. In figure 2, the ratio of the natural frequency w to 
the natural frequency wo obtained from elementary beam 
theory is shown as a function of the plan-form aspect ratio 
Lb for cross-sectional aspect ratios of 1.0, 3.6, and œ. 
The contribution of each of the secondary effects to the total 
reduction in the natural frequency for the cross-sectional 


aspect ratios 23.6 and 1.0 can be seen in figures 3 and 4, 


respectively. The cross-sectional aspect ratio of t- oœ 


corresponds to the limiting case of a beam where the effects 
of transverse shear deformation and longitudinal inertia are 
negligible and therefore the reduction in natural frequency 
is due entirely to shear lag. 

The dashed lines in figures 3 and 4 show the reduction in 
frequency due to the inclusion of the effect of only transverse 
shear deformation as obtained from reference 4. 

The long- and short-dash lines are calculated from the 
frequency equation (41) with krr=0 and consequently 


(a) First symmetrical mode. 


(b) Second symmetrical mode. 
Yiaurp 2.—Change in the natural frequency of a symmetrically vibrating free-free cylinder due to the inclusion of secondary effects. 


(c) Third symmetrical mode. 


t 


i zd symmetrical 


Effects included: 
—— — ——— Transverse shear (ref. 4) 
—— - ———— Transverse shear and shear lag 
——____—————- Transverse shear, shear lag, and 
2 longitudinal inertia 


5 6 10 I4 


L/b 


FiacurE 3.—Contribution of transverse shear deformation, shear lag, 
and longitudinal inertia to the reduction in natural frequency for 


=3.6. í 


represent the reduction in natural frequency when both shear 
lag and transverse shear deformation are taken into account. 
Thus the hatched area between the dashed and the long- and 
short-dash lines may be considered as showing the additional 
reduction in natural frequency when the influence of shear 
lag is considered. Finally, the solid lines are calculated with 
krr taken into account, and consequently the shaded area 
shows the additional ‘influence of longitudinal inertia in 
reducing the frequency. 


Examination of figures 3 and 4 and the curves for to œ in 
figure 2 shows that the influence of shear lag increases as the 
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1.0 
Mode: erry 
riid ree 
aranera eanont sateen 
8 
6 
wo 
A ES" 3d symmetrical 
t 
Effects included: 
—— —— —— Tronsverse sheor (ref. 4) 
——_-- Transverse shear and shear fag 
ò Transverse shear, shear lag, and 


longitudinal inertio 


O5 6 To) j4 
L/B 


FIGURE 4.—Contribution of transverse shear deformation, shear lag, 
and longitudinal inertia to the reduction in natural frequency for 


2=10. 
a 


cross-sectional aspect ratio increases; whereas the influence 
of transverse shear and longitudinal inertia decreases with 
increasing cross-sectional aspect ratio. Indeed, it appears 


from the results for t- 3.6 that for this aspect ratio the 


effects of ieagicddinal inertia may already be considered 
practically negligible. 

A word of caution concerning the interpretation of figures 
3 and 4 may be in order. Since in some cases the depth of 
the hatching increases with increasing L/b, it might appear, at 


_ first glance, that the shear-lag effect increases with increasing 


plan-form aspect ratio. However, if the additional effects 
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of shear lag are considered on a percentage basis with the 
dashed line as a base, it will be found that shear-lag effects 
actually reduce in percentage with increasing 0/b. A similar 
criterion should be used in judging the influence of longi- 
tudinal inertia. 


CONCLUDING REMARKS 


The numerical calculations show that secondary effects 
have appreciable influence on the natural frequencies of 
rectangular box beams of uniform wall thickness. These 
results constitute an indication of the probable inadequacy of 
elementary beam theory for the vibration analysis of actual 
aircraft structures of the monocoque and semimonocoque 
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type and emphasize the need for practical calculation pro- 
cedures for such structures that would take into account 
transverse shear deformation, shear lag, and, when necessary, 
longitudinal inertia. The general solutions presented for 
cylinders of uniform thickness, as well as the numerical results 
for rectangular box beams, should be useful in the assessment 
of the accuracy of any procedure of this kind that may be 
developed. 


Y 


LANGLÐY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
Lanceray FreLD, VA., January 21, 1962. 


APPENDIX 
TRANSFORMATION OF PARAMETERS 


Expressions for As/A, I, and K?.—Ifsin 6 is expanded into 
a Fourier series 
2nNT8 


(Al) 


the Fourier coefficients A, are the same as those defined in 
equation (24); that is, 


A,==Q sin 6 di OE (A2) 
Pp P 


The effective shear area Ag (eq. (25)) can now be written 
as a function of the Fourier series expansion for sin 0 as 


As=hie( È 
ahs 


With the use of the appropriate orthogonality conditions, 
equation (A3) becomes, after the integration is performed, 


A 

A, 2 
ee 
Ag 1 = 
TOR ee A? -o 
A 2 n=1,2,3 


2 
» sin “228 ds (A3) 


As= pt > Ap = 
Basas 
or 


(A4) 
The moment of inertia I of a cylinder is defined as (see 
fig. 1) 


io f ” ft da= AT? (A5) 


where ¥ is the y-distance to the center of gravity of the cross 
section and is given by 


(A6) 
But 
(A7) 


or 
(A8) 


(A9) 


* Inr 


With the use of equations (A8) and (A9), the expression 
for J in equation (A5) becomes 


_Ap’ As 
8 Sa n=]1,2,3 n? 


(A10) 


With the series expansion for J in equation (A10), the param- 
eter K?, as defined in equation (29), becomes 


| ee 
wut eS A: 


| m As akths nt a) 


Transformation of expression for N;—In equation (37) 
N; was defined as 


Nee oS. a_i 
New ge toes © Aoroa aB ER 2 


The infinite series that appears in this expression converges 
as A,” and therefore is a relatively slowly converging series. 
In order to increase its-rate of convergence, the following 
transformations are made.: 

By adding. and subtracting A,?/K* inside the infinite 
summation in equation (A12) and using equation (A4), the 
equation simplifies to 


O PPA pa & A,? ie | 
N= 16 Ag P neas kg Beta 9 he (A18) 
By adding and subtracting A,?/K?n? inside the infinite 


summation in equation (A13) and using equation (A11), the 
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expression for N; can be transformed to 


E Pr A 2 A, has 
Ni=gq BiT a, Bike’ 25 Res BFR 2? (A14) 


The infinite series in equation (A14) converges as A,?/n* and therefore is considerably more quickly convergent than the 
series in equations (A12) and (A13), which converge as A,’ and A,?/n’, respectively. 

Closed form of NW, for cylindrical beams of rectangular cross section.—For a cylindrical beam of rectangular cross 
section, with dimensions as shown in figure 2, it is possible to write the expression for N; in a closed form. The param- 
= eters for such a cross section become 


Agsg=4at ’ 
A=4(a+ b)i=pt 
A,=0 (n even) (A15) 
oe cos 2nrb (n odd) 
nr p 


With equations (A15) the parameter N; shown in equation (A12) becomes 


2 2NnT O ; 

vr È p < pil, ` 
Mi ShF TEE a naded papp 2 eae) 

K? 
or 
cog ene 
”2 2 -œ o 

AE PEE ST E p AREE N E = (M17 


N= 
i 8k Bka? a ibs papm pate Oe BAr? 


Each of the infinite summations in equation (A17) can now be written in closed form as shown in reference 5, and the 
closed expression for N; then becomes z 


ke 
sinh = £3 B z) 
Ba K p aK wa L., 
Ne=ges a2 meee + tanh 258 B, |)—= ks (A18) 
2K? 
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